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Monotone  Scmiflows  Generated  by  Functional  Differential  Equations 

Hal  Smith* 

Lefschetz  Center  for  Dynamical  Systems 
Brown  University 
Providence,  Rhode  Island  02912 

Recently,  there  has  been  a  considerable  advance  in  our  understanding  of  the 
qualitative  as  well  as  the  asymptotic  behavior  of  semiflows  on  partially  ordered 
spaces  which  preserve  the  partial  ordering.  In  large  part,  this  advance  is  due  to 
recent  work  of  MAV.  Hirsch  [9-12]  and  H.  Matano  [14-16].  The  most  striking  result, 
due  to  Hirsch  [10,11,12],  implies  that  "almost  every"  precompact  orbit  converges  to 
the  set  of  equilibria  (under  suitable  hypotheses,  see  section  4).  These  results  have 
been  applied  to  ordinary  differential  equations  in  Rn  (see  e.g.  [10,23])  where  the 
well-known  Kamke  theorem  applies  and  to  nonlinear  reaction  diffusion  systems 
with  quasimonotone  reaction  term  (see  e.g.  [11,12,14,15,23])  where  maximum 
principles  apply.  The  aim  of  this  paper  is  to  develop  the  machinery  necessary  to 
apply  the  above  mentioned  results  to  the  class  of  functional  differential  equations, 
FDE’s,  and  to  investigate  the  qualitative  behavior  of  the  subclass  of  these 
equations  which  generate  an  order  preserving  semiflow. 

More  precisely,  we  consider  the  FDE 

(0.1)  x’(t)  -  f(xt)  ,  •  -  d- 

at 

where  f:  C([-r,0],Rn)  -•  Rn  and  xt  denotes  the  element  of  C  =  C([-r,0],Rn),  the  space 
of  continuous  maps  of  [-r,0]  into  Rn,  given  by  x{(0)  *  x(t+8),  -r  <  0  *  0.  For  more 
details  concerning  FDE’s  we  refer  the  reader  to  the  text  [6].  Assume  that  (0.1), 
together  with  the  initial  data  Xq  ■  0  €  C,  gives  rise  to  a  unique  solution  on  [0,o), 
a  >  0,  which  we  denote  by  x(t,$)  or  xt($)  depending  on  whether  we  view  the 
solution  in  Rn  or  C.  Then,  under  suitable  conditions,  the  collection  of  maps, 
*On  leave  from  Arizona  State  University. 


$  -•  xt($),  is  a  local  semiflow  on  C.  If  4>,4>  €  C,  we  write  4>  $  4>  (0  <  0)  if  the 
indicated  inequality  holds  pointwise,  with  the  usual  (componentwise)  partial 
ordering  on  Rn.  The  semiflow  is  order  preserving  (we  will  say  that  f  is 
cooperative)  if  whenever  4>  S  4>  we  have  xt(4>)  «  xt(4>)  for  all  t,  0  <  t  <  minfo^o^}. 
Sufficient  conditions  for  f  to  be  cooperative  appear  not  to  be  well-known.  After, 
obtaining  such  conditions,  this  author  found  references  to  work  of  K.  Kunisch 
and  W.  Schappacher  [25],  R.H.  Martin  [13]  and  Y.  Ohta  [18]  who  had  earlier 
obtained  the  same  sufficient  conditions.  Most  likely,  others  before  them  have 
obtained  the  following  sufficient  condition 

(H)  Whenever  $  «  4>  and  ^(0)  =  4lj(0)  it  follows  that  fj(4>)  «  fj(40- 

For  those  familiar  with  the  Kamke  (quasimonotone)  condition  for  ordinary 

differential  equations,  (H)  will  seem  quite  natural,  it  reduces  to  the  Kamke 

condition. 

The  order  preserving  property  of  a  semiflow  is  not  sufficient  for  the 
strong  result  of  Hirsch  mentioned  above;  one  requires  strongly  order  preserving 
semiflows,  that  is,  if  4>  <  0,  $  *  4»  then  xt($)  <  xt(40,  at  least  for  all  large  t  (it  is 
only  reasonable  to  expect  such  an  inequality  for  t  >  r). 

In  section  2  of  this  paper  we  develop  sufficient  conditions  for  (0.1)  to 
generate  a  strongly  order  preserving  semiflow  (we  say,  in  this  case,  that  f  is 
cooperative  and  irreducible). 

In  section  3,  we  consider  the  stability  of  steady  states  of  cooperative 
FDE’s  and  the  existence  of  connecting  orbits  between  steady  states. 

In  section  4,  we  state  the  relevant  results  of  Hirsch  which  apply  in  our 

setting. 

The  main  result  of  our  work  can  be  roughly  summarized  as  follows.  Let 
f  be  cooperative  and  irreducible  (see  section  2)  and  assume  all  orbits  of  (0.1)  are 


precompact  (e.g.  all  orbits  are  bounded  and  f  maps  bounded  sets  to  bounded  sets). 
Then,  for  a  dense  set  of  initial  conditions  for  (0.1),  the  qualitative  behavior  of  the 
solutions  of  (0.1)  is  the  same  as  for  the  ordinary  differential  equation 


(0.2)  x’(t)  =  F(x(t)) 

F(x)  =  f(x) 


where  x  -*  x  is  the  inclusion  of  Rn  into  C  given  by  x(8)  =  x.  Equations  (0.1)  and 
(0.2)  have  the  same  steady  states  and  the  stability  properties  of  a  steady  state  x  of 
(0.1)  are  the  same  as  for  the  steady  state  x  of  (0.2).  Moreover,  Equation  (0.2)  is 
cooperative  in  the  sense  of  Hirsch  [9]  and  hence  there  are  simple  tests  for  stability 
and  instability  of  steady  states  of  (0.2)  (see  [23]  and  Corollary  3.2).  These  rather 
striking  results  represent  a  considerable  improvement  in  the  connection  made 
between  (0.1)  and  (0.2)  by  R.H.  Martin  [13]. 

The  results  of  the  first  four  sections  of  this  paper  are  applied  in  section 
5  of  this  paper  to  a  model  of  biochemical  feedback  in  protein  synthesis  which  goes 
back  to  Goodwin  [4]  and  has  been  the  object  of  much  study  [13,20,21,24] 
particularly  in  the  nondelay  case.  Under  very  reasonable  hypotheses,  we  obtain  an 
essentially  complete  picture  of  the  qualitative  behavior  of  the  solutions  of 
the  model  equations. 

It  is  convenient  to  establish  some  notation  here.  Let  R"  be  the  cone  of 
non-negative  vectors  in  Rn.  If  x,y  €  Rn  we  write  x  <  y(x  <  y)  if  xj  <  yj(xj  <  yj) 
for  1  <  i  <  n.  Let  (ej,...,en)  denote  the  standard  basis  in  Rn  and  let  l=ej-«2+...+en. 
If  x  <  y,  we  write  [x,y]  *  (z€Rn:  x  <  z  <  y).  Set  N  ■  (l,2,...,n).  Let  r  *  0  and 
C  *  C([-r,0],Rn)  be  the  Banach  space  of  continuous  functions  mappings  $:[-r,0]-*Rn 
with  supremum  norm.  If  $,tp  €  C,  we  write  0  <  4>(4>  <  <W  in  case  the  indicated 
inequality  holds  at  each  point  of  [-r,0].  If  4>  <  ip,  we  write  [$>,<[>]  =  (BeC:  <p  «  8  <  ip). 


V 


Let  C+  =  {0€C:  $  >  0}.  Let  denote  the  inclusion  Rn  -  C([-r,0],Rn)  by  x  -»  x, 

A 

x(8)  =  x,  6^-r,0],  Denote  the  space  of  functions  of  bounded  variation  on  [-r,0]  by 
BV[-r,0].  If  A  and  B  are  subsets  ot  a  linear  space,  then  A+B  =  {a+b:  aeA.beB}  and 
similarly  for  A-B.  If  X  and  Y  are  Banach  spaces,  let  L(X,Y)  denote  the  space  of 
bounded  linear  maps  from  X  to  Y.  Let  A  and  B  be  nxn  matrices.  We  write  A«B 
if  the  equality  holds  componentwise.  The  matrix  A  is  irreducible  if  it  does  not 
leave  invariant  any  proper  subspace  of  Rn  spanned  by  a  subset  of  the  standard 
basis  elements.  We  write  s(A)  =  max  Re  X  where  X  runs  over  the  eigenvalues  of  A, 
the  stability  modulus  of  A.  If  a  €  Rn  we  write  diag(a)  for  the  nxn  diagonal 
matrix  with  aj  in  the  (i,i)  entry. 


Let  ft  be  an  open  subset  of  R  x  C  and  f:  ft  -»  Rn  be  continuous.  Consider 

the  FDE 

(1.1)  x’(t)  =  f(t,xt). 

We  assume  throughout  this  paper  that  solutions  of  the  initial  value  problem  (1.1) 
together  with  xtQ  =  <J>,  for  (tg,tf>)  €  ft,  are  unique.  If  f  is  Lipschitz  continuous  in  <p 
on  compact  subsets  of  ft  then  uniqueness  holds  [6],  We  write  x(t,tg,<f>,f)(xt(tg,0,f)) 
for  the  solution  of  the  initial  value  problem  and  we  drop  the  f  when  no  confusion 
results. 

The  results  of  this  section  have  probably  been  proved  by  many  authors. 
However  we  are  only  aware  of  the  work  of  contained  in  [13], [18],  and  [25].  These 
authors  proved  both  results  of  this  section. 

Consider  the  hypothesis: 

(H):  If  (t,</>),(t,<W  €  n ,  <t>  i  4>  and  tf>j(0)  -  (/^(O)  for  some  i,  then 

fi(t,*)  S  fj(t,0). 

The  main  result  of  this  section  is  the  following. 

Proposition  1.1:  Let  f,g  :  n  -*  Rn  be  as  above  and  assume  either  f  or  g  satisfies 
(H).  Assume  f(t,<f>)  «  g(t,$)  for  all  (t,$)  €  a  If  (tg,<f>),(tg,<[»)  €  ft  with  <p  S  <p  then 
x(t,t0,<|>,f)  «  x(t,tQ,i^,g) 
for  all  t  >  tg  for  which  both  are  defined. 

Proof:  Assume  that  f  satisfies  (H).  Let  g€(t,$>)  =  g(t,4>)  +  el  and  4>e  =  <l>+€  1  for 
e  >  0.  If  x(t,tg,0,g)  is  defined  on  [tg-r,tj]  for  some  tj  >  tg  then  x(t,tg,i[>€,g£)  is 
defined  on  [tg-r,tj]  for  all  sufficiently  small  positive  e  by  Theorem  2.2  of  [6].  We 
will  show  that  x(t,tg,$,f)  <  x(t,tg,4>€ ,g€ )  on  [tg-r,tj]  for  small  positive  e.  The 
proposition  will  then  follow  by  letting  e  -*  0  and  applying  Theorem  2.2  of  [6]. 
Suppose  the  above  assertion  is  false.  Then  there  exists  a  small  positive  e  for 


which  x(t,tg,t|i€,g€)  is  defined  on  [tg-r,tj]and  a  t’€(tg,tj]  such  that 

x(t,t0,&f)  <  x(t,t0,i^€,g£)  on  [t0-r,t’)  and  Xj(t\tg,&f)  =  Xj(t\to,0€,g£)  for  some  value 

of  i.  Clearly  xj(t\tg,&f)  *  X;(t\tg,i/>£,g£).  But 

X-(t\t0,t|>£,g£)  =  gi(t’,Xt,(t0,4>€,g€»  +  6 
*  fi(t’,Xt.(t0,^,ge))  +  € 

>  fi(t’,Xt»(tO,0€,g€)) 

>  fjft’.Xj.ftg^.f)) 

=  X|(t’,tO,0,f ) 

where  the  latter  inequality  follows  from  (H).  This  contradiction  implies  that  such 
a  t’  cannot  exist  and  establishes  the  above  assertion  (and  the  proposition). 

Suppose  that  ft  =  R  x  U  where  U  is  an  open  subset  of  C  containing  C+. 
The  following  invariance  result  will  be  useful.  The  proof  involves  ideas  similar  to 
those  used  in  the  proof  of  Proposition  1.1  and  is  therefore  omitted. 

Proposition  1.2:  Assume  that  whenever  <P  e  C+  with  <^(0)  =  0  and  t  e  R,  fj(t,0)^O. 
If  <P  €  C+  and  tg  €  R  then  x(t,tg,4>)  >  0  for  all  t  >  tg  in  the  maximal  interval  of 
existence. 

Proposition  1.2  is  a  very  special  case  of  much  more  general  invariance 
result  due  to  G.  Seifert  [19]. 


The  ultimate  goal  of  this  section  is  to  find  sufficient  conditions  for  the 
autonomous  FDE 

(2.1)  x’(t)  =  f(xt),  f  :  C  -  Rn 

to  have  the  property  that  whenever  0,0  €  C  are  distinct  with  0  (  0, 

then  xt(O,0)  S  xt(0 ,0)  for  t  >  0 

and  xt(O,0)  <  xt(O,0)  for  t  >  tQ  >  0, 

t0  sufficiently  large  (independent  of  0  and  0).  The  following  example  will  show 
us  that  we  need  to  modify  slightly  our  notion  of  the  state  space  of  (2.1). 

Consider  the  initial  value  problem 
x  j  (t)  =  axj(t)  +  bx2(t-l/2) 

t  *  0 

x2(t)  =  cxj(t-l)  +  dx2(t) 

Xj(0)  HO 

,  -1  <  e  <  o 

x2(0)  =  02(6) 

supp  <P2  C  t-1,-2/3],  <t>2  *  0,  02  p  0. 

The  initial  value  problem  can  easily  be  integrated  by  steps  and  one  finds 
x(t)  =  0  for  t  >  -2/3.  The  problem  is  that  x2(t-l/2),  t  >  0,  "never  sees"  the  support 
of  02  so  x’j(t)  =  0.  Hence,  although  the  initial  condition  0  =  (O,02)  >  0  satisfies 
0  0,  we  nevertheless  have  x(t,O,0)  agreeing  with  the  identically  zero  solution. 

The  source  of  the  failure  of  the  solution  operator  to  distinguish  between 
the  two  initial  conditions  0  and  0  is  our  (implicit)  choice  of  the  state  space  as 
C  »  C([-1,0],R.2).  We  will  show  that  this  pathology  can  be  removed  by  taking  our 
state  space  to  be  C^jjy2j  s  C([-1,0],R)  *  C([-1/2,0],R).  To  this  end,  consider  the 
above  linear  equation  where  b  >  0,  c  >  0.  Let  (0j,02)  €  C^j  j/2)  and  assume 
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(0j,$2)  *  0.  We  can  associate  to  (#j,$2)  =  0  a  solution  x(t,O,0)  and  it  is  not 
difficult  to  see  that  x(t,0,$)  >  0  for  t  >  3/2  (the  worst  case  is  where  the  support  of 
$1  belongs  to,  say,  [-€,-«/2]  for  small  e  and  $2  =  0)-  *n  addition,  for  each  small 
€  >  0,  one  can  select  a  nontrivial  $  as  above  such  that  xj(t,0,$)  =  0  on  [0,3/2-€]  but 
x(t,0,$)  >  0  for  t  ?  3/2. 

It  is  worth  emphasizing  this  last  point.  Namely  that,  with  our  new 
choice  of  state  space,  if  <p  >  0  and  4>  *  0,  the  solution  x(t,0,$)  has  the  property  that 
x(t,0,$)  >  0  for  t  £  3/2.  In  other  words,  it  takes  time  for  nontrivial  nonnegative 
initial  conditions  to  become  a  positive  state.  Fortunately,  we  can  bound  from 
above  this  time  lapse  independently  of  the  initial  conditions. 

Motivated  by  this  example,  we  develop  some  notation.  Let  r  €  Rj  , 
r  =  (rj,r2,...,rn),  |r|  =  max  r;  and  let 

Cr  =  C([-r  j,0],R)  x  C([-r2,0],R)  x  ...  x  C([-rn,0],R) 

C  =  C([-|r|,0],Rn). 

Given  an  element  <p  =  ($j,...,$n)  €  Cf  we  will  identify  4>  with  an  element  of  C, 
which  again  we  label  4>,  by  the  correspondence 

#8)  =  (*i(e)-..,*n(0)}  where  :  [-|r|,0] -R 

$j(9)  -rj  «  6  S  0 

is  defined  by  ?j(0)  = 

- 1  r  |  t  8  <  -rj 

This  inclusion,  I,  of  Cf  into  C  is  injective  but  not  surjective.  Similarly,  given  an 
element  0  of  C  we  will  identify  $  with  an  element  of  Cr,  which  we  label  4>  again, 
by  4>  =  (0j,...,0n)  €  Cr  ,  <t>i  -  ft),  -ej  . 

[*ri>0] 

This  latter  map,  R,  of  C  into  Cr  is  surjective  but  of  course,  not  injective.  We 
write  C*  for  the  cone  of  n-tuples  of  nonnegative  functions  and  if  4>,(p  €  C*  we 
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write  <p  (  ip  (0  <  0)  if  0-0  €  C*  (0-0  €  interior  C+).  Note  that  we  use  the  same 
symbol  for  inequality  in  both  C  and  Cr  This  should  not  cause  a  problem  since 
0  (  0  in  Cr  if  and  only  if  0  *  0  in  C  (upon  making  identifications  as  above). 
Finally,  let  Ker(R)  =  (0€C  :  0;(0)  -  0,  >rj  <  6  <  0,  1  «  i  <  n}. 

It  will  be  convenient  to  define  the  semigroup  on  C  generated  by  the 
trivial  FDE  :  x’=0.  For  T  £  0,  define  S(T)  :  C  -•  C  by 

'  #6)  - 1  r  |  S  0  S  -T 

[S(T)«(0)  =  if  T  <  |r| 

<p( 0)  -T  (  0  s  0 

and 

[S(T)0]  =  0(0)  if  T  ?  |  r  |  . 


With  this  notation,  we  begin  by  considering  a  nonautonomous  linear 
equation.  Our  aim  is  to  establish  sufficient  conditions  for  certain  nonnegative 
nonzero  initial  conditions  to  give  rise  to  eventually  positive  solutions. 

Consider  the  linear  FDE 
(2.2)  x’(t)  =  L(t,xt) 

where  L(t, •)  :  C  -  Rn  is  a  bounded  linear  map  and  t  -*  L(t,-)  is  continuous  from  1R 
to  L(C,Rn),  the  space  of  bounded  linear  maps  of  C  into  Rn.  The  following 
hypotheses  will  be  required 


In  order  to  get  another  perspective  on  the  hypotheses  (K)  and  (I), 
consider  the  standard  representation  of  L  as  L  =  (Lj,...,Ln)  where 
n  0 


M(t,« 


=  E  J 

j=1  -|r 


0j(8) 


1  i  i  <  n 


in  which  n- j  :  RxR  -»  R  satisfies 


nij(6,t)  =  njj(0,t)  0  )  0 

hij(9,t)  =  j(- 1 r | ,t)  SO  9<-|r| 

r>i j(-,t)  €  BV[-|r|,0] 

njj(  -  ,t)  is  continuous  from  the  left  on  (-|r|,0). 


If  we  set 

a j(t)  s  nH(0,t)  -  nH(0-,t),  1  «  i  <  n 


hjiO.t) 


^(6,0, 


0  <  0 
0)0 


then 


0  0 

Li(t,0)  =  ajOWjto)  +  J  ^(ejdQOjjce.t)  +  e  J  ^j<e)dei7ij(e,t) 

i*i  ,|r| 


-  r 


=  ai(t)  0j(O)  +  Lj(t,«  . 


It  is  easy  to  see  that  (K.)  holds  if  and  only  if  i?jj  and  >7jj,  i*j,  are  nondecreasing  in 

0  for  fixed  t.  Hence  (K)  holds  if  and  only  if  L(t, •)  :  C+  -*  R”  . 

It  is  not  difficult  to  see  that  due  to  the  continuity  of  the  map  t  -»  L(t, ■), 


a j(t)  is  continuous  in  t. 


Finally,  we  note  that  the  matrix  A(L)(t)  has  entries 


0 

(2.3)  A(L)(t)ij  =  J  dQh^e.t)  *  njjCO.t)  . 

*|r| 

The  next  several  results  concern  solutions  x(t,tQ,0)  of  (2.2).  We  begin  by 
showing  that  (K)  implies  that  if  a  solution  x(t,tQ,0)  (0  €  C+)  ever  has  a  positive 
component  then  ever  after  that  component  is  positive. 

Lemma  2.1:  Let  (K)  hold,  0  €  C+  and  x(t)  =  x(t,tQ,0),  t  *  tQ,  satisfy  (2.2). 

If  Xj(tj)  >  0  for  some  tj  >  tQ  then  xj(t)  >  0  for  t  >  tj. 

Proof:  Note  that  (K)  implies  the  hypotheses  of  Proposition  1.2  so  x(t,tQ,0)  >  0  for 

t  *  tQ.  Now 

Xj(t)  ■  -aj(t)xj(t)  +  Lj(t,xt)  >  -a j(t)x j(t). 

Hence,  if  xj(tj)  >  0  then  xj(t)  >  0  for  t  >  tj  by  a  standard  differential  inequality 
argument. 

The  next  lemma  is  the  rationale  for  assuming  (R).  It  guarantees 
"ignition”  of  some  component  of  x(t,tQ,0). 

Lemma  2.2:  Let  (R)  hold.  If  0  €  C+,  0  f  Ker(R),  tQ  €  R  then  there  exists 

i  €  N  and  t  €  [tQ,tQ+|r|)  such  that  Xj(t,tQ,0)  >  0. 

Proof:  If  0(0)  *  0  then  the  conclusion  follows  by  taking  i  €  N  such  that 

0j(O)  >  0  and  t  =  tQ. 


Assume  #0)  *  0.  By  (R),  there  exists  i  €  N  and  T  €  [0,  |r|)  such  that 
Lj(t0+T,S(T)*)  >  0.  Since  [S(T)fl(0)  =  0,  Li(t0+T,S(T)$)  =  Li(t0+T,S(T)*).  Clearly 
xt(t 0,<t>)  >  S(t-t0)<#>  for  t  >  t0  so  Lj(to+T,xto+T(to,0))  >  Lj(t0+T,S(T)$)  >  0.  Hence 
xtQ+T(tQ^)  *  0  and  this  implies  our  result. 

Putting  together  (K.)  and  (R),  we  see  that  if  4>  €  C+,  <p  f  Ker(R)  then 
some  component  of  x(t,tg,$)  becomes  and  remains  positive.  To  "turn  on"  the  other 
components  is  the  job  of  (I). 

Proposition  2.3:  Let  (K),(I)  and  (R)  hold.  If  Q  €  C+,  4>  f  Ker(R),  tQ  €  R  then 
x(t,to,0)  >  0  for  t  >  tQ+n|r|  . 

Proof:  By  Lemmas  2.1  and  2.2,  there  exists  an  i  €  N  such  that  xj(t)  >  0  for 

t  *  tQ+|r|.  It  follows  that  there  exists  a(t)  >  0  for  t  >  t0+2  | r  |  such  that  xt*a(t)ej 
for  t  >  t0+2 1  r  | .  Hence  L(t,xt)  >  a(t)L(t,ej)  for  t  *  t0+2 1  r  | .  Now  since 
A(L)(t)  =  -[diag  a(t)]  +  A(L)(t),  A(L)(t)  is  irreducible.  It  follows  that  there  exists 
j*i  such  that  Lj(tQ+2 |r |,ej)  >  0.  Either  Xj(tQ+2 1 r | )  >  0  or  Xj(tg+2|r|)  =  0  and 

x’j(t0+2|r|)  -  Lj(t0+2|r|,xto+2|r|) 

=  Lj(t0+2|r|,xto+2r)  >  «(t0+2  |r  |)  Lj(t0+2  |r  |,Cj)  >  0. 

But  xj  (tQ+2  |  r  | )  =  0  and  Xj(tQ-t-2 1  r  | )  >  0  contradicts  Xj  (t)  *  0  for  t  >  tQ.  Hence  we 
must  have  Xj(tg+2|r|  >  0  and  by  Lemma  2.2,  Xj(t)  >  0  for  t  >  tg+2|r|. 

It  follows  that  there  exists  8(t)  >  0  for  t  >  tg+3|r|  such  that  xt  >  a(t)ej+B(t)ej 
for  t  >  tg+3 1 r  | .  Hence  L(t,xt)  >  ^tJLCt.ej)  +  B(t)L(t,ej)  =  A(L)(t)(a(t)ei+0(t)ej). 
Now  A(L)(tg+3  |r  |)  is  irreducible  so  A(L)(tg+3  |r  |)(a(tQ+3  | r  |)ej+0(to+3  | r  |)ej)  ek>0 
for  some  k*{i,j}.  The  reasoning  is  as  follows.  Since  A(L)(tj)  (tj  ■  tg+3|r|)  is 
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irreducible  A(L)(t3>  does  not  leave  invariant  the  span  of  ej  and  ej.  It  follows 
that  A(L)(t3)(Moc(t3)ej  +  X0(t-j)  ej)ek  >  0  for  some  choice  of  jt€(0,l),  V€(0,1)  and 
k€N-(i,j).  Since  A(L)(t3)  is  nonnegative  the  assertion  above  follows.  Now,  if 
xk(t3)  *  0  then 

xic(l3)  =  Ut3,xt3)  >  A(L)(t3)(a(t3)ei+e(t3)ej)  •  eR  >  0. 

But  this  is  incompatible  with  xk(t)  *  0  for  t  >  tQ,  hence  xk(t3)  >  0.  By  Lemma 
2.2,  xk(t)  >  0  for  t  >  tQ+3|r|. 

Continuing  in  this  manner,  we  obtain  x(t)  >  0  for  t  £  t0+n|r|. 

Recall  our  earlier  example  where 
L(t,0)  =  col(a  0j  (O)+b02(- 1  /2),c0j  (- 1  )+d02(O)) 
where  b  >  0  and  c  >  0.  Then  (K)  holds  with  L(t,0)  =  col(b02(-l/2),c0j(-l)).  The 
matrix  A(L)(t)  is  given  by 

a  b 

A(L)(t)  = 

c  d 

b 

and  is  clearly  irreducible  since  both  b  and  c  are  positive.  It  is  easy  to  check  that 
(R)  holds  for  this  example  if  r  =  (1,1/2)  determines  the  restriction  map  R.  Notice 
that  (R)  fails  if  "0  f  Ker(R)"  is  omitted  from  (R). 

We  now  turn  to  the  main  goal  of  this  section.  Let  r  €  Rj  and 
C  -  C([-|r |,0],Rn).  Consider  (2.1)  where  f  :  U  -  Rn  is  a  continuously 
differentiable  map  on  the  open  subset  U  of  C.  Assume  U  is  order  convex,  that  is, 
if  0,0  €  U  with  0  <  0  then  t0  +  (l-t)0  €  U  for  0  <  t  <  1. 
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Definition:  f  is  cooperative  in  U  if  for  every  <p  €  U,  L  =  df(40  satisfies  (K). 

f  is  cooperative  and  irreducible  in  U  if  f  is  cooperative  in  U  and  the  following 
hold: 

(1)  If  4>  €  U  and  ip  €  C  with  $-ip  €  Ker(R)  then  ip  €  U  and 

m  =  m. 

(2)  For  all  0  €  U,  L  =  df (0)  satisifies  (I)  and  (R)  with  the 
provision  that  i  and  T  are  independent  of  ip. 

We  will  write  x(t,$)  (x{($))  for  x(t,0,$>)  €  Rn  (xt(0,<>)  €  C).  Thus  x(t»  is 
the  solution  of  (2.1)  satisfying  \q(4>)  =  4>.  The  next  lemma,  together  with 
Proposition  1.1  implies  that  if  f  is  cooperative  in  U  and  0  (  <p  then  xt(4>)  i  xt(ip) 
for  t  >  0  on  their  common  interval  of  existence. 

Lemma  2.4:  If  f  is  cooperative  in  U  then  f  satisfies  (H)  of  section  1. 

Proof:  Suppose  4>,ip  €  U,  4>  4  <p  and  ^j(O)  =  ipj(O).  Since  U  is  order  convex  and  f 

is  continuously  differentiable 

1 

=  Jdfj(s<fH-(l-s)<W(<M)ds  *  0 
0 

where  the  inequality  holds  since  the  integrand  is  pointwise  nonnegative  by  (K). 

The  next  result  is  the  main  theorem  of  this  section.  It  will  allow  us  to 
define  a  local  semi-flow  on  Cr  which  is  eventually  strongly  monotone. 

Theorem  2.5:  Let  f  be  cooperative  and  irreducible  in  U.  If  and  ip  are 

elements  of  U  with  $  4  </>,  «{*-$  f  Ker(R)  and  [0,o),  0  <  o  <  «  is  the  intersection  of 
the  maximal  intervals  of  existence  of  x(t,$)  and  x(t,<M,  then 
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x(t,0)  <  x(t ,0)  for  n  |  r  |  <  t  <  a  . 

Proof:  By  Lemma  2.4  and  Proposition  1.1  we  have  x(t ,0)  *  x(t,0)  on  0  «  t  <  a  . 

Now 

1 

x(t,40  -  x(t,0)  =  J  d0x(t,s0+(l-s)0)(0-0)ds. 

0 

If  5  €  U  and  0  €  C,  d^xft.OB  =  y(t,0)  satisfies  the  linear  variational  equation 
y’(t)  =  df(xt(5))yt  ,  y0  *  0 

(See  Theorem  4.1  of  [6].)  Let  L(t, •)  =  df(x{(0).  It  is  apparent  that  L  satisfies  (K) 
and  (I).  In  order  to  see  that  (R)  holds  for  L,  suppose  ft  €  C,  n  f  Ker(R)  and 
H(0)  =  0.  By  our  assumption  that  f  is  cooperative  and  irreducible,  there  exists  i 
and  T  )  0  such  that  df(Oi(S(r)jx)  >  0  for  all  [  6U.  In  particular,  if  j;  »  xTU), 
?€U,  we  have  Lj(t,S(t)m)  >  0  which  is  just  (R). 

We  can  now  apply  Proposition  2.3:  if  0  €  C+,  0  f  Ker(R)  then 
y(t,0)  =  d^x(t,5)B  >  0  for  t  >  n|r|.  Now  <!>-$  €  C+  and  0-0  f  Ker(R)  so  for  each 
fixed  s,  the  integrand  above  is  positive  for  t  J  n|r|.  Hence  the  integral  above  is 
positive  and  we  have  proved  the  theorem. 

If  f  is  cooperative  and  irreducible  in  U,  we  would  like  to  think  of  the 
state  space  for  (2.1)  as  Cr  rather  than  C.  The  reason  for  this  is  simple.  In  order 
to  apply  a  very  powerful  result  of  Hirsch  for  strongly  monotone  flows  [11,12]  (see 
section  4)  we  need  the  semi-flow  generated  by  the  maps  0  -*  xt(0)  to  have  the 
strong  monotonicity  property  that  if  0  <  0  and  0*0  then  xt(0)  <  xt(0),  at  least 
for  large  t.  Theorem  2.5  does  not  give  us  this  in  C.  However,  in  Cr,  we  have  this 
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property  as  we  shall  show.  First  note  that  if  i p,4i  €  C  and  €  Ker(R)  then 
x(t,$)  =  x(t,ip)  for  t  >  0.  Here  is  where  (1)  of  the  definition  of  cooperative  and 
irreducible  plays  a  role  (a  proof  can  be  given  following  the  proof  of  Theorem  2.5). 
Consider  the  collection  of  maps 

*t 

<p  €  Cr  - >  R  xt(I<f>)  ,  t  Z  0 

defined  on  open  subsets  of  U,  depending  on  t.  It  is  apparent  that  {4>t}  defines  a 
local  semi-flow  in  the  sense  of  Hirsch  [11,12],  Moreover,  if  <p,<p  €  Cr,  <P  $  ip,  <p  *  <l>, 
then  clearly  14>  <  I<l>  and  I4»-1 ^  Ker(R).  Hence,  by  Theorem  2.5,  xt(I</>)  >  x{(10) 
for  all  t  *  (n+l)|r|  for  which  both  exist  and  thus  <  4>t(0)  for  t  £  (n+l)|r|. 
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3.  Stability  of  steady  states  and  connecting  orbits. 

In  this  section  we  are  concerned  with  the  stability  of  a  steady  state  of 

the  FDE 

(3.1)  x’(t)  =  f(xt) 

where  f  is  a  continuously  differentiable  cooperative  map  f  :  U  -*  Rn,  U  an  open 
order  convex  subset  of  C,  r  e  Rj.  Suppose  there  exists  v  €  Rn  with  v  €  U  and 
f(v)  =  0 

and  consider  the  linear  variational  equation  about  the  steady  state  v: 

(3.2)  y’(t)  =  L(yt)  ,  L  =  df(C). 

Let  L  be  represented  as 


L  4>  =  (L  L„0) 
n  0 

Lj*  -  Z  I  ♦:(e)dnij(0) , 
j-‘  ~ I r | 


1  S  i  S  n 


njj  €  BV[-|r|,0]  ,  Ojj(.|r|)  =  0  ,  *..(%)  continuous  from  the  left  on 

(•|r|,0). 


The  characteristic  values  associated  with  the  linearized  equation  (3.2)  are  roots  of 
det  A(  X)  =  0 

(3.3)  A(X)  =  VI  -  A(V) 


A(V) 


U 


J  eX0dbij(8) 


-  r 


Define  the  stability  modulus  of  L  as 

s(L)  ■  max(Re  V  :  det  A(V)  =  0)  . 

The  stability  modulus  is  a  well-defined  quantity  since  for  any  B  e  R  there  are 
at  most  a  finite  number  of  zeros  of  det  A(V)  with  Re  V  >  3.  It  is  well  known  that 
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the  steady  state  v  is  asymptotically  stable  if  s(L)  <  0  and  is  unstable  if  s(L)  >  0. 

Since  f  is  cooperative  L  satisfies 
(K)  Whenever  <p  €  C+  and  <f>j(0)  =  0  ,  <p  *  0. 

Hence  the  bjj  satisfy 


j(®)  >  ‘  *  L  *s  nondecreasing  in  9  €  [-|r  ,0] 


(3.4)  nii(6)  = 


^(9)  ,  e  €  [-  r  ,0) 


^ii(O-)  ,9  =  0 


n- -( 0)  is  nondecreasing  in  9  €  [-  r  ,0]. 


Our  first  result  says  that  the  stability  of  a  steady  state  of  a  cooperative 
system  is  determined  by  a  real  "most  unstable"  characteristic  root 


Theorem  3.1:  Let  f  be  cooperative.  Then  s(L)  is  a  root  of  det  A(x)  =  0.  If  X 

is  a  characteristic  root  different  from  s(L)  then  Re  X  <  s(L).  If  L  satisfies  the 
hypotheses  of  Proposition  2.3  then  s(L)  is  a  simple  root  of  det  A(X)  =  0. 

Proof:  Consider  A(X)  for  real  values  of  X.  For  i*j,  it  follows  from  (3.3),(3.4) 

that  A(X)jj  ?  0.  Hence  A(X)  is  a  matrix  with  nonnegative  off  diagonal  elements. 
Similarly,  if  X^  i  X2,  then  A(X2)"  <  A(Xj)|j.  This  is  immediate  from  (3.3)  and 

9  ,  a  _ 

(3.4)  for  i*j.  For  i=j,  note  A(X)jj  =  hjjtO-n^O-)  +  j  exodnii(9)  and  hjj  is  non- 

•lrl 

decreasing  so  A(X2)-  <  A(Xj)jj.  Thus  we  have  that  Xj  «  X2  implies  A(X2)  »  A(Xj). 

Consider  what  happens  as  X  -*  ®.  From  (3.3)  we  have  that 

lim  Ajj(X)  =  ni i(0)-»?i |(0-),  hence  limA(X)  =  A(®)  =  (b: :(0)-n::(0-))  and  observe 
X-®  X-*  j  j 


A(®)jj  >0  i*j. 


V  <V 
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Let  us  write  s(A(X))  for  the  stability  modulus  of  the  matrix  A(X).  Since 
for  real  X,  A(X)  has  nonnegative  off  diagonal  elements,  it  is  known  that  s(A(X))  is 
an  eigenvalue  of  A(X).  Furthermore,  if  Xj  i  X2,  A^)  i  A(Xj)  and  so 
s(A(X2))  S  s(A(X j)).  It  is  also  known  that  s(  )  is  continuous,  thus  the  map 
X  -*  s(A(X))  is  a  nonincreasing  continuous  map  from  R  into  itself  which  has  a 
finite  limit  at  +®.  It  follows  that  there  exists  a  unique  value  of  X,  Xq,  for  which 
s(A(X0))  =  Xq.  We  will  show  that  s(L)  =  Xq. 

Since  s(A(Xq))  =  Xq  and  s(A(Xq))  is  an  eigenvalue  of  A(Xq),  it  follows  that 
det(s(A(X0))I  -  A(Xq))  =  det(X0I  -  A(XQ))  =  0 
and  hence  Xq  is  a  characteristic  root  of  (3.2).  Assume  X  is  a  real  characteristic 
root  of  (3.2)  so  that  det(XI-A(X))  =  0.  Then  X  is  an  eigenvalue  of  A(X)  so  XSs(A(X)). 
It  follows  that  X  <  Xq.  Indeed,  if  X  >  Xq  then  X  <  s(A(X))  S  s(A(Xq))  =  Xq  <  X. 
Thus  Xq  is  greater  than  or  equal  to  any  real  characteristic  root  of  (3.2). 

Let  T(t)0  =  yt(0)  denote  the  solution  of  (3.2)  with  <f>  €  C.  Then 
T(t)(C+)  C  C+  and  T(t)  is  a  compact  linear  operator  for  t  »  |r|.  That  is,  {T(t))t^Q 
is  a  strongly  continuous  positive  semigroup  consisting  of  compact  operators  for 
t  *  |  r  | .  Now  if  p(T(t))  denotes  the  spectral  radius  then  it  is  known  that  [6] 
p(T(t))  =  et,z,  u  =  max{Re  X  :  det  A(X)  =  0} . 

On  the  other  hand  the  set  of  characteristic  roots  is  precisely  the  spectrum  of  the 


infinitesimal  generator  of  (T(t))  (see  [6]).  It  has  been  shown  that  ti  belongs  to  the 
spectrum  of  the  generator  of  a  strongly  continuous  positive  semigroup,  [5].  In 
particular,  n  is  a  characteristic  root  so  n  =  Xq. 

Moreover,  by  Theorem  7.2  in  [8,  see  also  17],  Xq  is  the  only  characteristic 
root  X  satisfying  Re  X  ■  Xq. 


eigenvalue  of  T(t)  and  thus  Xq  is  a  simple  characteristic  root  [1],  This  completes 
our  proof  of  Theorem  3.1. 

The  proof  of  Theorem  3.1  yields  much  more  than  has  been  stated.  In  the 
next  few  remarks,  we  bring  out  other  consequences  of  the  proof. 

Remark  1:  There  exists  u  ?  0  in  Rn  such  that  y(t)  =  ues^)1  satisfies  the 

variational  equation  (3.2).  This  is  an  immediate  consequence  of  the  well  known 
fact  [2]  that  corresponding  to  s(A(Xq)(=s(L)=Xq)  there  is  a  nonnegative  eigenvector 
in  Rn.  If,  in  addition,  L  satisfies  the  hypotheses  of  Proposition  2.3,  then  one  can 
take  u  >  0.  Again,  this  follows  from  the  fact  that  A(Xq)  is  irreducible  which  in 
turn  follows  from  (I)  of  the  previous  section. 

Remark  2:  The  importance  of  the  fact  that  s(L)  is  itself  a  characteristic 

root  of  (3.2)  will  be  clear  to  anyone  who  has  had  experience  in  computing 
characteristic  roots  of  FDE’s  --  even  one  dimensional  equations.  They  are 
notoriously  difficult  to  find.  It  is  very  important  to  be  able  to  determine  stability 
by  only  considering  the  real  characteristic  roots  of  (3.2). 

Remark  3:  The  second  assertion  of  Theorem  3.1  is  especially  important  in 

the  context  of  bifurcation  theory.  One  naturally  asks  the  question  "how  can  a  one 
parameter  family  of  steady  state  of  a  parametrized  family  of  FDE’s  lose  stability 
at  a  critical  value  of  the  parameter?”  According  to  Theorem  3.1,  if  the  family 
consists  of  cooperative  maps,  the  answer  is  that  a  real  eigenvalue  must  change  sign 
giving  rise  to  a  steady  state  bifurcation  (generically)  and  an  exchange  of  stability. 
In  particular,  although  a  Hopf  bifurcation  is  not  precluded  for  cooperative  FDE’s 
(see  [21]),  a  steady  state  can  never  lose  stability  to  a  Hopf  bifurcation.  Put 
another  way,  a  local  Hopf  bifurcating  periodic  solution  cannot  be  asymptotically 
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stable  (see  also  Theorem  4.2  of  section  4). 


The  following  result,  a  corollary  to  the  proof  of  Theorem  3.1,  provides  a 
relatively  simple  test  for  stability  and  instability. 


Corollary  3.2:  s(L)  <  0  (s(L)  >  0)  if  and  only  if  s(A(0))  <  0  (s(A(0))  >  0). 

Moreover  s(A(0))  <  0  if  and  only  if 


(-1)J 


where  A(0)jj  =  >7^(0). 


A(0)j  j  .  .  .  A(0), j 
A(0)j!  .  .  .  A(0)jj  J 


>  0  ,  j  =  1,2 . n. 


Proof:  The  second  equivalence  is  well  known  [2].  Recall  s(L)  =  XQ,  the  unique 

root  of  X-s(A(X))  «  0.  But  X  -*  X-s(A(X))  is  an  increasing  continuous  function. 
Hence  XQ  <  0  if  and  only  if  0-s(A(0))  >  0,  that  is,  if  and  only  if  s(A(0))  <  0. 
Similarly  XQ  >  0  if  and  only  if  0-s(A(0))  <  0. 


Corollary  3.2  has  an  interest  interpretation.  It  is  saying  that  the  zero 
solution  of  the  linear  FDE  is  asymptotically  stable  or  unstable  according  as  the 
zero  solution  of  the  linear  ODE 
Z’  =  A(0)Z 

is  asymptotically  stable  or  unstable.  Notice  how  A(0)  is  obtained:  A(0)jj  =  hjj(0); 
the  magnitude  of  the  delays  are  completely  ignored!  There  is  an  appealing  way  to 


look  at  A(0),  that  is 
A(0)v  -  L(v) 


in  other  words,  A(0)  is  the  restriction  of  L  to  the  constant  functions.  This 
observation  leads  to  another  interesting  observation.  Consider  the  nonlinear 
ordinary  differential  equation 
(3.5)  x’(t)  =  F(x(t))  ,  F(x)  «  f(i)  . 

It  is  easy  to  see  that  (3.5)  is  cooperative  in  the  sense  of  Hirsch  [9]  (see  also  [23]). 
The  steady  states  of  (3.5)  and  (3.1)  are  identical.  Even  more,  the  stability  type  of 
a  steady  state  is  the  same  for  (3.5)  as  for  (3.1)!  One  only  need  check  that 

A 

0X 

DF(x)  =  df(£)  ^ —  df(x)(e  j,...,en)  *  A(0). 


Hence  as  far  as  steady  states  and  their  stability  goes,  one  can  trade  in  a 
cooperative  FDE  for  a  cooperative  ODE.  The  first  person  to  see  a  connection 
between  (3.5)  and  (3.1)  appears  to  have  been  R.H.  Martin  [13].  He  did  not  make 
the  observation  that  steady  states  have  the  same  stability  type,  however. 

If  L  satisfies  the  hypotheses  of  Proposition  2.3  and  s(L)  >  0  then  by  our 
earlier  remark,  there  exists  a  solution  of  (3.2)  of  the  form  y(t)  =  ues^L^,  u  >  0.  In 
terms  of  our  state  space  Cp  this  solution  gives  rise  to  a  monotone  orbit  t  -»  yt(u), 
t  €  R,  of  (3.2)  where  u  €  C*  is  given  by 


(3.6) 


u  =(uj,...,un) 

uj(0)  =  u^1^0  ,  -rj  <  0  <  0  . 


The  orbit  connects  the  trivial  solution  of  (3.2)  to 


This  "most  unstable" 


manifold  for  the  variational  equation  should  have  a  counterpart  for  the  nonlinear 
equation  (3.1).  This  is  precisely  the  content  of  the  next  theorem. 
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Thcorcm  3.3:  Let  f  be  cooperative  and  irreducible.  Suppose  f(v)  =  0, 

A  A 

s  «  s(df(v))  >  0  and  suppose  df  is  Lipschitz  continuous  in  a  neighborhood  of  v. 
Suppose  v  +  C+  belongs  to  the  domain  of  f  and  f  is  bounded  on  bounded  subsets 
of  U.  Then  there  exists  a  unique  C*  function  y^O,®)  -*  v  +  C*  satisfying 

(1)  y(r)  =  v  +  Tu  +  o(T)  as  T  -*  0 
where  u  >  0  is  as  ir  (3.6). 

(2)  xt(y(T))  =  y(estT)  ,  t  *  0,  T  >  0. 

(3)  0  (  Tj  S  t2  implies  y(Tj)  <  y(T2). 

(4)  Either  (a)  lim  ||y(T) ||  =  •  or  (b)  lim  y(T)  =  w  where 

J-4QO 

w  €  Rn,  w  >  v  and,  f(w)  =  0  and  s(df(w))  «  0. 

(5)  If  (4)(a)  holds,  then  for  all  *  v  ,  0  *  v,||xt(4>)||  -*  «®  as 
t  tends  to  the  right  hand  limit  of  the  maximal  interval 
of  existence  of  xt($).  If  (4)(b)  holds,  then  for  all 

0,  v  <  4>  «  w,  4>  *  v,  xt($)  -*  w  as  t  -*  “  . 

Theorem  3.3  says  that  the  monotone  curve  T  =  {y(T):T  )  0}  is  a  heteroclinic  orbit  of 
(3.1)  connecting  the  unstable  steady  state  v  to  the  steady  state  w  (or  ®).  In 
addition,  the  steady  state  w,  if  not  asymptotically  stable  (s(df(w))  «  0),  at  least 
attracts  all  initial  conditions  $  different  from  v  with  v  i  4>  <  w.  When  w  in 
Theorem  3.3  exists,  we  will  paraphrase  Theorem  3.3  by  simply  saying  that  there 
exists  a  monotone  increasing  trajectory  connecting  v  to  w. 

Actually,  we  have  only  stated  half  of  the  story.  If  v-C+  belongs  to  the 
domain  of  f  then  one  can  find  a  monotone  decreasing  trajectory  connecting  v  to  p 
(or  •)  with  the  obvious  changes  in  (l)-(5)  above. 

The  assumption  that  v+C+  belongs  to  the  domain  of  f  can  be 
significantly  weakened.  For  example,  if  it  is  known  that  there  exists  a  steady 
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state  w  with  w  >  v  and  [v,w]  belongs  to  the  domain  of  f  then  one  can  show  that 
r  C  [v,w]  . 

The  proof  of  Theorem  3.3  is  very  similar  to  the  proof  of  Theorem  2.7  in 
[23]  and  uses  Theorems  1.1  and  2.1  in  [22].  Hence  we  will  not  give  the  proof  here. 
It  should  be  remarked  that  Matano  [15,16]  has  stated  a  similar  result,  although 

A  A 

requiring  a  second  steady  state  w  >  v  but  less  smoothness  than  we  require. 


In  this  section  we  state  some  very  powerful  results  due  to  M.W.  Hirsch 
[9- 1 2 j  for  monotone  semi-flows  which  have  direct  application  in  our  setting.  We 
will  not  state  these  results  in  there  greatest  generality.  Let  X  be  a  separable 
Banach  space  and  K  a  cone  in  X,  that  is,  K.  is  a  nonempty  closed  subset  of  X  with 
the  closure  properties  R+  K  C  K,  K+K  c  K.  and  K  n(-K.)  =  {0}.  We  assume  that  K 
has  nonempty  interior,  K°.  We  write  x  <  y  (x  <  y)  if  and  only  if  y-x  belongs  to 
K.  (K°).  Let  U  be  an  open  order  convex  subset  of  X.  Let  $>  =  {4>t}t be  a  local 
semi-flow  on  U  (see  e.g.  [11,16]),  in  particular,  each  $t  is  a  continuous  map  on  a 
subset  Ut  of  U  and  the  semigroup  property  *t*s  -  ®s+t  holds  under  appropriate 
conditions.  We  say  that  $  is  a  monotone  flow  if  each  <frt  is  monotone:  x,y  €  Ut 
and  x  <  y  implies  <X>t(x)  «  C*t(y);  ♦  is  strongly  monotone  if  each  4>t  is  strongly 
monotone:  x  <  y,  x  *  y,  t>0  implies  4>t(x)  <  4>t(y). 

As  an  example,  consider  a  cooperative  FDE  defined  on  an  order  convex 

*t 

subset  of  Cr,  r  €  R" .  Then  X  *  Cr  and  the  local  semi-flow  is  given  by  £__>xt(0). 
This  semi-flow  is  monotone  in  the  above  sense.  If  f  is  cooperative  and  irreducible 
then  4>t  is  eventually  strongly  monotone  but  not  strongly  monotone.  That  is,  there 
exists  T  >  0  such  that  if  x  <  y,  x  #  y  then  $t(x)  <  4>t(y)  for  t  >  T  (provided 
4>t(x),4>t(y)  exist).  For  FDE’s  we  may  take  T  =  (n+l)|r|  by  Theorem  2.5.  Note  that 
T  is  independent  of  x  or  y.  The  results  of  Hirsch  which  we  state  below  for 
strongly  monotone  flows  are  true  for  what  we  term  eventually  strongly  monotone 
flows  (see  also  Matano  [16]). 

The  first  two  results  require  only  a  monotone  flow  (See  [11],  Theorem 


2.3  and  Corollary  2.4.) 
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Theorem  4.1:  Let  x  €  U  be  such  that  i9+(x)  -  {4>t(x):  t  >  0}  is  compact  in  U. 

Assume  for  some  real  T  >  0,  ♦-j-(x)  >  x  (4>t(x)  <  x).  Then  4>t(x)  converges  to  an 
equilibrium  as  t  -*  *. 

Theorem  4.2:  A  monotone  flow  does  not  have  an  attracting  periodic  orbit. 

By  an  equilibrium,  we  mean  a  point  y  €  U  such  that  4>t(y)  *  y  for  all  t*0. 
By  periodic  orbit  we  mean  a  non-constant  closed  orbit.  Such  an  orbit  is  attracting 
if  it  attracts  an  open  set.  It  is  not  difficult  to  see  that  Theorem  4.1  implies 
Theorem  4.2. 

The  following  results  require  the  flow  to  be  strongly  monotone.  For 
simplicity,  we  assume  that  the  domain  of  4>t  is  U  for  every  t  >  0  and  that  all 
orbits  have  compact  closure  in  U.  The  following  result  is  a  special  case  of 
Theorems  5.2  and  5.5  in  [11]. 

Theorem  4,3:  Let  S  be  a  totally  ordered  subset  of  U  and  let  E  be  the  set  of 

equilibria.  Assume  E  consists  of  isolated  points.  Then  the  subset  of  points  of  S, 
the  orbits  of  which  do  not  converge  to  a  point  of  E,  is  at  most  countably  infinite. 

In  particular.  Theorem  4.3  implies  that  a  dense  set  of  points  of  U  have 
convergent  orbits. 

The  following  result  of  Hirsch  [12]  sharpens  the  conclusion  of  Theorem 
4.3  at  the  expense,  of  course,  of  additional  hypotheses.  It  does  not  require  that  X 
be  separable  though  (see  Theorem  10.1  in  [12]). 


Theorem  4.4:  Let  Xg  be  an  open  positively  invariant  set  for  {4>t)tjg  and 

suppose  Xg  contains  a  compact  attractor  which  attracts  points  of  Xg.  Assume  E  is 
a  finite  set.  Then  there  is  a  dense  open  subset  Dg  of  Xg  consisting  of 
"asymptotically  stable"  points  x  for  which  u(x)  C  E. 

By  u(x),  we  mean  that  Omega  limit  set,  n  cl(  u  4>t(x)),  of  the  orbit 

T>0  t>T 

through  x.  A  compact  attractor  K  is  a  compact  invariant  set  which  attracts  a 
neighborhood  of  itself.  In  Theorem  4.4,  we  assume  u(x)  C  K  for  all  x  €  Xg.  We 
remark  that  EOK  is  nonempty  (see  [11]). 

We  now  make  use  of  our  remarks  concluding  section  2,  together  with 
Theorem  4.4  and  Lemma  2.2  in  [6]. 

Theorem  4.5:  Assume  (2.1)  is  cooperative  and  irreducible  in  U  and  that  f 

maps  bounded  subsets  of  U  to  bounded  subsets  of  Rn.  Assume  that  U  is  positively 
invariant  for  (2.1),  there  is  a  closed  bounded  subset  B  of  U  such  that  for  all 
0  €  U,  u(<f>)  €  B  and  the  map  0  -»  xt($)  maps  bounded  sets  to  bounded  sets  for  each 
t  *  0.  Then  B  contains  a  compact  attractor  which  attracts  points  of  U.  If  E  is  a 
finite  set  consisting  of  nondegenerate  equilibria,  then  the  union  of  the  basins  of 
attraction  of  the  equilibria  e  with  s(df(e))  <  0  is  an  open  dense  subset  of  U. 

Proof:  By  Lemma  2.2  in  [6],  B  contains  a  compact  attractor  which  attracts  points 

of  U.  By  Theorem  4.4,  there  is  a  dense  open  subset  Ug  of  U  consisting  of 
asymptotically  stable  points  0  for  which  u($)  c  E.  But  if  0  is  an  asymptotically 
stable  point  with  u($)  -  (c),  e  €  E,  then  necessarily  s(df(e))  <  0.  Hence  Ug  is  a 
subset  of  the  union  of  basins  of  attraction  of  equilibria  e  with  s(df(e))  <  0. 


As  an  application  of  the  ideas  in  the  previous  section  we  consider  a 
mathematical  model  of  biochemical  feedback  in  protein  synthesis  developed  by 
Goodwin  [4]  and  which  has  been  the  object  of  much  study  [13,20,21,24],  Wc  refer 
the  interested  reader  to  [4]  for  details  concerning  the  model.  The  quantities 
xj,X2,...,xn  denote  concentrations  of  mRNA  (xj),  various  enzyme-intermediates 
(x2,...,xn_j)  and  a  final  product  protein  (xn)  in  a  sequence  of  first  order  reactions 
Xj  -*  x2  -*  ...  -*  xn.  The  product  protein  is  assumed  to  induce  (positive  feedback) 
the  transcription  of  mRNA. 

The  equations  derived  from  the  model  are  given  below  as  (5.1).  In  (5.1), 
we  have  used  the  notation  x;t  to  denote  the  function  x;t(0)  =  x:(t+0),  -r;  (  0  <  0. 

Jil  J  J 

x’j(t)  -  h(L  x  t)  -  ajXj(t) 

(5.1) 

x’j(t)  =  Lj.jxj.j  j  -  «jXj(t),  2  <  j  <  n 
0 

Lj4)  -  J  WB)dr)i(0) 

'ri 

’’i  •  [-rj,0]  -*R  is  nondecreasing 

^(-rj)  =  0,  TTj(O)  =  1,  Hj(0)  >0  for  0  >  -rj 

h  :  R+  -*  R+  is  a  smooth  function  with  h’(u)  >  0  for  u  >  0 
oj  >  0,  1  <  i  <  n  . 

R.H.  Martin  has  considered  equation  (5.1)  in  the  special  case  that  h(u)  =  u/l+u  and 
Lj0j  -  ^j(-rj)  in  [13].  We  have  borrowed  techniques  from  Martin  as  well  as  from 
previous  authors  (e.g.  Selgrade  [20],  Hirsch  [10])  who  have  treated  simpler  versions 
of  (5.1).  Allwright  [24]  also  treats  a  version  of  (5.1). 

A  brief  word  on  the  assumptions  listed  below  (5.1)  is  appropriate. 
Concerning  the  integrators  Pj,  we  set  hj(-rj)  -  0  merely  as  a  normalization  (recall 


t\ j  and  i?|+c,  c  €  R,  give  equivalent  integrals).  Our  requirement  that  hj(0)  =  1  can 
always  be  achieved  by  a  scaling  of  the  xj  and  an  appropriate  redefinition  of  h.  If 
our  requirement  that  hj(0)  >  0  for  6  >  -rj  does  not  hold  then  Oj(0)  =  0  on  some 
subinterval  [-rj,-Sj]  and  hence  one  could  replace  -rj  by  -Sj.  This  assumption,  given 
rjj(O)  =  l,  is  thus  without  loss  of  generality.  On  the  other  hand,  it  is  essential  for 
the  proper  choice  of  state  space,  namely  Cr,  r  =  (r  1^2.-, r„)  >  0,  for  (5.1)  to  be 
cooperative  and  irreducible.  We  require  h’  to  be  locally  Lipschitz  in  order  to 
apply  Theorem  3.3. 

We  write  the  right  hand  side  (5.1)  as  f(xt)  where 

=  col(h(Ln^n)-«1^1(0),L1^1-«2^2(0)! . V^n-rCnV0)) 

and  f  :  C  -*  Rn. 

Since  h  and  Lj  are  nondecreasing  in  their  arguments,  f  satisfies  (H),  in 
fact,  f  is  cooperative.  Note  that  this  would  not  be  the  case  if  -ajx(t-Tj)  replaced 
-aix(t).  In  addition  f  satisfies  the  hypotheses  of  Proposition  1.2  so  that  the  flow 
<P  -*  x{($),  t  >  0  leaves  C+  positively  invariant  (so  long  as  solutions  are  defined). 

Our  first  task  is  to  insure  that  solutions  with  nonnegative  initial 
conditions  are  globally  defined  and  to  seek  conditions  for  boundedness  of  these 
solutions.  The  following  inequality  will  be  useful. 

h(u) 

(5.2)  lim  sup -  =  a  <  *  . 

u  -•  +®  u 

Lemma  5.1:  Suppose  (5.2)  holds.  If  0  €  C+  then  x(t,$)  is  defined  for  t>0 

and  x(t,$)*0.  If 

(5.3)  a  <  fla- 

then  $+($)  -  (xt(0):  t  >  0)  has  compact  closure. 


Proof:  By  (5.2)  there  exists  b  >  0  such  that  h(u)  *  au+b  for  u  *  0.  Let 
g(u)  =  au+b  and  write  x(t,$,h),  x(t,$,g)  for  the  solution  of  (5.1)  and  (5.1)  with  g 
replacing  h,  respectively.  The  solution  x(t,$,g),  £  €  C+  is  globally  defined  since  it 
satisfies  a  linear  FDE.  By  Proposition  1.1  and  Proposition  1.2  we  have 
0  i  x(t,$,h)  <  x(t,0,g)  on  the  maximal  interval  of  existence  of  x(t,tf>,h).  Now  f 
maps  bounded  subsets  of  C+  to  bounded  subsets  of  Rn  so  by  Theorem  3.2  in  [6], 
x(t,$,h)  can  fail  to  be  globally  defined  only  by  becoming  unbounded,  a  possibility 
that  our  inequality  precludes.  It  follows  that  x(t,&h)  is  defined  for  t  *  0. 

Now  suppose  (5.3)  holds.  The  function  x(t,tf>,g)  satisfies  a  linear  non- 
homogeneous  equation  which  we  write 

(5.4)  Z’(t)  *  LZt  +  bej 

where  L$  -  col(aLn0n-a10j(O),L|^j-a2^2(°) . Ln-l^n-l'^^n^0^-  Wc  wiu  show  that 

s(L)  <  0  by  making  use  of  Corollary  3.2.  First  we  calculate  A(0)  (see  (3.3)) 

-ocj  0  .  .  .  0  a 

1  -02  0  .  .  .  0 

(5.5)  A(0)  =  0  1  -o3  0  .  .  0 

,  0  •  •  •  0  1  -«n 

Note  that  the  subdiagonal  of  ones  arises  due  to  our  normalization  nj(0)  *  1.  Now 
it  is  easy  to  check  that  the  n  principal  minors  of  Corollary  3.2  have  the  correct 
sign  if  (5.3)  holds.  Hence  the  zero  solution  of  the  linear  homogeneous  equation 
Z’(t)  *  LZt  is  asymptotically  stable.  But  this  implies  that  the  positive  steady  state 
xn  ■  b/(oj  O2...otn-a),  xj  »  ^j+l  xj+l>  j*lt-»n-l,  of  the  nonhomogeneous  equation  (5.4) 
is  globally  attracting.  The  boundedness  of  x(t,$)  *  x(t,$,h)  follows  from  the 


comparison  0  <  x(t,0,h)  S  x(t,0,g).  The  compactness  of  the  orbit  -Q +(0)  follows  from 
the  fact  that  it  is  bounded,  that  f  maps  bounded  sets  to  bounded  sets,  and  the 
Ascoli-Arzela  theorem. 

We  have  observed  that  f  is  cooperative.  In  fact,  it  is  cooperative  and 
irreducible  as  we  now  show.  For  0  €  C+, 

df(M  =  col(eLn^n-a101(O),L1^ra202(O),...,Ln.1«n.1-an0n(O)) 

0  =  8(0)  =  h’(Ln0n)  >  0. 

In  order  to  check  (I),  we  need  to  show  that  the  matrix  A(df(0))  is  irreducible  for 
each  0  €  C+.  But  A(df(0))  is  identical  to  the  matrix  A(0)  in  (5.5)  except  that 
a=6(0).  Hence  A(df(0))  is  irreducible  if  and  only  if  8(0)  ■  h’(Ln0n)  >  0. 

Consider  the  requirement  (R)  (tg=0).  Let  0  e  C+,  0  f  Ker  R  and  0(O)=O. 
Then  for  some  j  and  some  0g^-rj,O],  0j(0o)  >  0.  Set  i=j+l  where  we  agree  that 
n+l  =  l  and  set  T  *  rj+0g  >  0.  Note  that  both  i  and  T  depend  only  on  0  and  not  on 

0 

0.  If  j<n  then  df(0)i(S(T)0)  =  Lj(S(T)0)  =  J  (S(T)0)j(6)d»7j(0),  where 

'rj 

(S(T)0)j(-r j)  -  0j(0g)  >  0.  Since  we  assumed  Oj(8)  >  hj(-rj)  =  0  for  0  >  -rj,  it 
follows  that  the  integral  is  positive  since  the  support  of  (S(T)0)j  and  the  support 
of  dOj  overlap.  If  j=n  and  8(0)  >  0,  a  similar  calculation  shows  df(0)j(S(T)0)>O. 
We  have  proved  the  following 

Lemma  5.2:  f  is  cooperative  and  irreducible  in  C+. 

We  now  turn  to  the  question  of  existence  and  stability  of  steady  states  of 
(5.1).  We  assume  that  (5.3)  holds.  The  steady  states  of  (5.1),  x,  x  €  Rj  are 


solutions  of 
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$ 


h(xn)  =  «jxi 

xj-l  "  “jxj  >  2  «  J  «  n- 

and  are  in  one-to-one  correspondence  with  solutions  of 
h(xn)  -  «!a2  ...  =  0  ,  x„  >  0  . 

We  assume  that 

(5.6)  0  is  a  regular  value  of  q(u)  s  h(u)  -  aj  ...  anu  . 

It  follows  from  (5.3)  and  (5.6)  that  there  exists  at  least  one  steady  state  and  there 
are  finitely  many  steady  states.  Moreover,  the  steady  states  are  totally  ordered 


a  i  a Am 

X*  <  <  ...  <  X*** 


The  stability  of  x*  depends,  according  to  Corollary  3.2,  on  the  stability  of  the 
matrix 


-«1 


0  1 


The  principal  minors  alternate  in  sign  as  in  Corollary  3.2  if  and  only  if 
h’(x{,)  <  aj...an.  The  reverse  inequality  implies  s(Aj)  >  0  as  is  easy  to  check.  We 
have  proved  part  of  the  following 


1  -<*2  0 


0  h’(x|j) 

0 
0 


1 


Theorem  5.3:  Assume  (5.3)  and  (5.6)  hold.  Then  xm,xm*2,xm*4,...  are 

asymptotically  stable  and  xm'1,xm'^,...  are  unstable.  There  exists  a  monotone 
increasing  orbit  connecting  xm'*  to  xm  and  a  monotone  decreasing  orbit 


connecting  to  xm**  to  xm'2.  An  identical  assertion  holds  for  the  other  unstable 


steady  states  xm*^,.., 


Proof:  (5.3)  and  (5.6)  imply  that  h’(xJ^)-aj...on  <  0,  thus  xm  is  asymptotically 

stable.  By  (5.6)  the  sign  of  h’(x|1)-a|...an  must  alternate  with  i.  The  existence  of 
connecting  orbits  follows  from  Theorem  3.3. 


The  main  result  of  this  section  follows.  We  write  B(x‘ )  for  the  domain 
of  attraction  of  x‘,  i=m,m-2,m-4,...  .  That  is,  B(x')  =  {(|)€C+:  u($>)  =  {x*}}. 


Theorem  5.4:  Assume  (5.3)  and  (5.6)  hold.  Then  [0,xm]  attracts  all  orbits  of 

(5.1)  with  4>€C+.  If  m=l  then  B(x*)  =  C+.  If  m  >  1,  then  u  B(xm'^)  is  open  and 

j*0 


dense  in  C+. 


Proof:  In  Lemma  5.1  we  showed  that  every  orbit  has  compact  closure  in  C+ 

and  in  fact  is  attracted  to  a  bounded  subset  of  C+,  namely,  to  [0,y  ]  where 
y€Rn  is  the  equilibrium  of  the  linear  inhomogeneous  comparison  equation. 
Hence  the  hypotheses  of  Theorem  4.5  are  satisfied.  We  only  need  to  verify 
the  conclusion  B(x^)  =  C+  when  m=l.  We  first  observe  that  x^  +  C+  belongs 
to  B(xJ).  In  fact,  if  <t>  Z  x1  then  there  exists  <l>  with  x1  S  <t>  <  such  that 
i/>  €  B(x')  (use  Theorem  4.4  or  Theorem  4.5).  But  then  €  B(xJ)  since 

x1  «  xt(tf>)  S  xt(ip)  for  t  i  0.  Similarly,  if  <P  €  C+  and  ^Sx1  then  <f>  €  B(x*). 
Hence,  if  <p  €  C+  there  exist  $j,$2  €  Btx1)  with  $>j  (  M  ^  and  hence  <P  € 
B(x').  This  completes  our  proof. 


One  might  conjecture  that  every  orbit  of  (5.1)  with  nonnegative  initial 
condition  converges  to  an  x1  even  when  m  >  1.  However,  this  is,  in  general, 
false.  There  can  be  periodic  orbits  of  (5.3),  necessarily  unstable  (Theorem  4.2). 
Indeed,  Selgrade  [21]  has  shown  for  the  ODE  version  (L:xit  *  x;(t))  of  (5.1) 

J  Jtl  J 
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that  a  Hopf  bifurcation  can  occur  at  an  unstable  steady  state. 

It  is  worth  mentioning  that  by  Theorem  4.5,  (5.1)  possesses  a  compact 
attractor  in  [0,xm]  C  C+.  This  compact  attractor,  J,  is  invariant  and  connected. 

A  1 

Of  course,  when  m=l,  J  =  {x  },  but  if  m  >  1  then  J  consists  of  the  equilibria,  their 
connecting  orbits  described  in  Theorem  5.3,  and  the  unstable  manifolds  of  the 
equilibria  xm'*,xm'^,...  .  Any  exotic,  but  necessarily  unstable,  dynamics  together 
with  its  attracting  set  must  be  connected  in  J  by  the  unstable  manifolds  of  the 
unstable  equilibria. 
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